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FLUTTER OF ELASTICALLY SUPPORTED ORTHOTROPIC 
PANELS INCLUDING THE EFFECTS OF FLOW ANGLE 


By James Wayne Sawyer 
Langley Research Center 

SUMMARY 

A theoretical panel flutter analysis and computer program was developed which 
is capable of analyzing an orthotropic panel with in-plane loads, at various angles of 
cross flow and with two opposite edges mounted on flexible supports and the other two 
edges simply supported. The resulting modal analysis uses linear piston-theory aero- 
dynamics and includes both aerodynamic and structural damping. Calculations were made 
for typical panels with no in-plane forces showing the effect on flutter of various deflec- 
tional edge- support conditions, flow angle, and aerodynamic and structural damping. The 
results show that large reductions in dynamic pressure for flutter are possible with only 
small changes in flow angle, and the reduction is greater for elastically supported panels 
than for simply supported panels. Aerodynamic damping has a significant stabilizing 
effect at all flow angles except zero, whereas structural damping has negligible effect. 

INTRODUCTION 

Interest in reentry and sustained high-speed flight has resulted in the need for a 
system to protect the vehicle's primary structure from high surface temperatures. Sev- 
eral types of thermal protection systems are being developed, including both metallic and 
nonmetallic heat shields. (See, for instance, refs. 1 to 4.) One type of thermal protection 
system of interest is shown in figure 1 and consists of a flexible thermal insulation blanket 
covered by a wavy metallic surface (corrugated or corrugation stiffened) which is attached 
to the primary structure by flexible supports. The supports must be flexible enough to 
allow thermal expansion without imposing large stresses but strong enough to transmit 
airloads to the primary structure. The wavy surface geometry required to accommodate 
thermal strain results in a panel with highly orthotropic stiffnesses which in service will 
experience flow over its surface at angles other than parallel to the corrugations. As a 
result of the flexible supports and flow angle, panel flutter becomes an important factor 
to consider in the design of such a panel. This dependence on flow angle and flexible sup- 
ports has been shown experimentally in reference 5 in which orthotropic panels mounted 
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on flexible supports were shown to experience large reductions in critical flutter dynamic 
pressures for only small changes in flow angle. 

Previous theoretical investigators have studied the influence of boundary conditions 
and flow angle on flutter of orthotropic panels separately. For example, references 2 
and 6 to 9 have studied the effects of flexible supports on flutter for panels at flow angles 
of 0° and 90°, and references 10 to 13 have studied the effects of arbitrary flow angles on 
flutter for panels with classical (simply supported or clamped) boundary conditions. No 
theoretical panel flutter studies have been made, however, for panels mounted on flexible 
supports and at arbitrary flow angles. 

The present investigation was undertaken to develop an analysis and computer pro- 
gram that is capable of analyzing an orthotropic panel at various flow angles and with 
two opposite edges mounted on flexible supports and the other two edges simply supported. 
The resulting modal analysis uses linear piston-theory aerodynamics and includes in-plane 
loads and both aerodynamic and structural damping. The aerodynamic theory treats the 
surface as a flat plate and does not take into account wavy surfaces such as shown in 
figure 1. The number of modes in each direction can be varied to insure that converged 
solutions are obtained. Results are presented for typical panels with no in-plane loads 
showing the effect on flutter of various edge- support conditions, flow angle, and aerody- 
namic and structural damping. Flutter mode shapes and frequencies are presented for a 
typical panel at various flow angles. 


SYMBOLS 

A Fourier series coefficients (see eq. (6)) 

A quantity defined by equation (15a) 

a panel length (x-direction) 

B frequency parameter defined by equation (15b) 

b panel width (y-direction) 

C mn constant of integration 

c free-stream speed of sound 

Dmn constant of integration 
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panel bending stiffness in x- direction 


D, 


xy 


D l = 


panel bending stiffness in y-direction 


panel twisting stiffness 




°2 = 


1 


Mx^y 


d 12 ~ ^xy + ^x^2 

g a aerodynamic damping coefficient (see eqs. (13)) 

bending structural damping coefficient 
membrane structural damping coefficient 


Sb 

g 


m 


K d >K r jK t deflectional, rotational, and torsional spring constants, respectively, per 
unit length 

K d ,K r ,K(. nondimen s ional deflectional, rotational, and torsional spring constants, 
respectively 


M 


m 


N x ,N y 


n 


r,s 


Mach number 

number of half-waves in streamwise direction 

uniform in-plane loads per unit length in x- and y-directions, respectively 

number of half-waves in cross-stream direction 

dynamic pressure of airstream 

integers 

time 



w 


lateral deflection of panel 


w 


mn 


X 


x,y 

x,y 

a 

tt mn 

J3 

y 

e } 6 

A 

X 

x cr 

Px>Py 


V 

P 

* 

U) 


lateral deflection of panel describing shape of natural mode of vibration 

function describing shape of natural mode of vibration in x-direction 
(see eqs. (16) and (17)) 

Cartesian coordinates of panel (see fig. 2) 

fixed coordinates based on stream direction (see fig. 2) 

complex frequency coefficient 

complex natural frequency coefficient 

compressibility factor 

panel mass per unit area 

parameters associated with characteristic roots (see eqs. (18) and (19)) 
local cross-flow angle, deg 
dynamic -pres sure parameter, 2qa^/^Dj 
dynamic-pressure parameter at flutter 

Poisson’s ratio associated with curvature in y- and x-directions, respectively 
nondimensional coordinate, x/a 
nondimensional coordinate, y/b 
free-stream air density 
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real part of complex-frequency coefficient 
imaginary part of complex-frequency coefficient 



to r 


fundamental frequency of simply supported beam, 
radians per second 



Subscripts: 

m, r number of half-waves in x-direction 

n, s number of half-waves in y-direction 

ANALYSIS 


Flutter Equation 

A flutter analysis is presented for a flat orthotropic panel with supersonic flow at 
Mach number, M, over one surface at arbitrary cross-flow angle, A. The panel and 
coordinate system are shown in figure 2. The panel edges at x = ± are supported by 
deflectional, rotational, and torsional springs, and the other two edges (y = 0 and b) are 
simply supported. The panel is loaded by uniform in-plane loads, N x and Ny, which 
are considered positive in compression. The aerodynamic loading is given by linear 
piston theory (which assumes a smooth exterior surface) and includes aerodynamic 
damping. The governing differential equation assuming linear small-deflection theory is 
(see ref. 10) 


( 


1 +ig h )[D 1 + 2D 12 - + D 2 

' b '\ 1 9x4 12 ax20y2 2 gy4 

9t P 0X 


V 


“ill j 


dx* 


oy 


(1 


where 

D x 

D 1 

1 “ Mx^y 

D 

b 2 = — X_ 
1 “ Mx^y 

= D xy + 
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and g m are bending and membrane structural damping coefficients (see ref. 14), 
and pc is the aerodynamic damping term. After transforming the aerodynamic 
loading terms to panel coordinates (see fig. 2) and nondimensionalizing, equation (1) 
becomes 



The boundary conditions at the simply supported edges are 

w(x,0,t) = w(x,l,t) = (x,0,t) = (x,l,t) = 0 (3) 

dT) z 

At the spring -supported edges, the boundary conditions obtained from reference 15 are 



(4) 


(5) 
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and K d> K r , and K t are the respective deflectional, rotational, and torsional spring 
constants per unit length. 


Flutter Solution 

Assume the solution of equation (2) to be in the form 


w = 


oO OO 


Z Z ^mn w mn e 
m=l n=l 


at 


( 6 ) 


where a = xf/+ ico. The natural mode shapes, w mn , satisfy the boundary conditions, 
equations (3), (4), and (5), and the following natural vibration equation obtained by setting 
the dynamic terms in equation (2) equal to zero: 


(‘ 


+ 






+ w 
^ U u mn w mn 


.filler 


mn w mn = 0 


(7) 


The complex frequency coefficient, a nin , is the damped natural frequency coefficient of 
the mnth mode. 

Substituting the expression for w given by equation (6) into equation (2) and 
making use of equation (7) results in the following: 



If w mn is taken in the form 

Wjjyj = X mn (£ )(sin n7TTj) (9) 

the term sin m r\ satisfies the boundary conditions, equation (3), directly. A Galerkin 
type solution is then obtained by substituting the expression for w mn given by equation (9) 
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into equation (8), multiplying through by X rs sin snrj, and integrating over the total area. 
Performing the integration and rearranging the terms results in the following set of 
simultaneous equations: 



+ 



mrns 

mrns 


l 

7T 4 b S =1 



sin A = 0 


(n / s) 


( 10 ) 


where 


1/2 aX^ 

^mrns “ f- 1/2 X rs dx 

_ ,1/2 

P mrns " 1/2 X mn X rs dx 


and 



7T 4 Dj 
y a 4 








( 11 ) 

( 12 ) 


(13) 


For a nontrivial solution to equation (10), the determinant of the coefficients of A mn 
must equal zero. A standard eigenvalue routine for a square complex matrix was used 
to calculate the eigenvalues of equation (10) which render the determinant equal to zero. 
The complex frequency coefficients, a , are the eigenvalues. Flutter is considered to 
occur for that value of the dynamic-pressure parameter, X, for which one of the imagi- 
nary frequencies changes sign (see ref. 14). 


Natural Modes and Frequencies 

The natural (no-flow) mode shapes, X mn , and frequency coefficients, a mn , are 
needed to determine the coefficients in equation (10). These are obtained from the 
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differential equation (7) and boundary conditions, equations (4) and (5). Substituting equa- 
tion (9) into equation (7) results in the following differential equation, in 


XU + Att 2 x" - Bv 4 X m „ = 0 
mn mn mn 


(14) 


where 



^. 2 (g) 2 ^2 n 2 

Djjr 2 \W 

[ gm| N y a4 /anV» D 2 ya V p^a 

is b/ Djb 2 ?!- 2 ' D l Dj^l + ig b ) Djjr 4 ^ + ig b ) 


(15a) 


(15b) 


Appropriate mode shapes which satisfy the boundary conditions, equations (4) and (5), are 
given as follows (see ref. 7): 

For even modes 

^mn ~ ^mn cos * 1 e mn^ + ®mn cos (16) 

For odd modes 

x mn _ 6 mn sinh + l^mn s ^ n ®mn^ (17) 

The relationship between e mn and 6 mn and the coefficients of equation (14) are 
obtained upon substitution of equation (16) or equation (17) into equation (14): 
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In order to satisfy the boundary conditions, the expressions for X mn given by 
equations (16) and (17) are substituted into equation (9) which is then substituted into 
equations (4) and (5). Solution to the resulting two simultaneous equations provides a 
transcendental equation in e mn and 5 mn . 

Method of Computation 

In order to obtain the aerodynamic coefficients of the flutter determinant (eq. (10)), 
mode shapes and frequencies are required. These are obtained for a given value of A 
by choosing B and solving for e mn and 6^ in equations (18) and (19). These values 
are then substituted into the transcendental equation. If the transcendental equation is 
not satisfied, B is varied and the process is repeated. With the final values of e mn 
and 6,^, the mode shapes, X mr) , are defined within an arbitrary constant and the fre- 
quency coefficients, a mn , are obtained from equation (15b). Since the coefficients 
as defined by equations (15) are complex, the resulting frequencies are complex 
quantities . 

When flow angle is not a variable, general flutter calculations can be made that 
cover a wide range of panels. Two geometry -stiffness parameters that are usually spec- 
ified in the literature (see ref. 6) and supply all the information needed to make general 
calculations for zero flow angle are 



These parameters are sufficiently general that one solution is applicable to panels with 
various length-width and stiffness ratios. With the inclusion of flow angle, however, 
equation (10) shows that the panel length-width ratio, a/b, must also be specified. By 
specifying a/b, the parameters A^. and C become functions only of the ratios of 
flexural stiffness properties, D 12 /Dj and D 12 /d 2 - Thus, flutter calculations must be 
made for the specific panel of interest. Flutter results are presented in the next section 
for panels with a wide range of stiffness parameters. 

RESULTS AND DISCUSSION 
Range of Parameters 

Calculations were made for a wide range of flow angularity, boundary conditions, 
and panel characteristics. Flow angle was varied between 0° and 90°, and calculations 
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were made for deflectional spring stiffness, K^, of 10, 100, 500, and In order to 
keep the variables to a reasonable number, however, in-plane loads, N x and N y , and 
spring supports, K r and K{., were set equal to zero. Aerodynamic and structural 
damping coefficients, g a , g b , and g m , were set equal to zero for all panels except for 
one representative panel, which was used to study the effects of damping on flutter. Cal- 
culations were made for of 0.05, 0.25, and 0,50; *> 12/^2 °* 20 > 50, and 500; 

and (a/b) 2 of 1.0 and 10. This range of parameters was chosen because they include 
most orthotropic panels that have been investigated either experimentally or theoretically 
in the literature (see, for instance, ref. 16) and most panels that are of current interest. 

Convergence 

In normal-mode flutter analyses, one of the questions that naturally arises is how 
many modal terms are required to obtain a converged solution. Previous studies have 
shown that, for flow parallel to a panel edge (zero flow angle), one modal term normal to 
that edge is sufficient in flutter analyses, and the number of modal terms required par- 
allel to the edge depends on the degree of orthotrophy (see, for example, ref. 7). When 
the flow angle is not zero, however, several modal terms may be required in both direc- 
tions, but the number is dependent on panel orientation and boundary conditions. Using 
more modal terms than are necessary results in a considerable waste in computer run 
time and storage and may make the solution prohibitively costly. 

Convergence data are given in figures 3(a), 3(b), and 3(c) where X^/3 i s shown as 
a function of the number of modal terms for a typical orthotropic panel with various 
deflectional spring edge support conditions and at flow angles of 0°, 90°, and 10°, respec- 
tively. Calculations are made for D 12^1 = ^*5, D 12,/ D 2 = (a/b) 2 = 1.0. For 

the extreme orientations (A = 0° and 90°), calculations indicate only one modal term is 
required in the direction perpendicular to the flow to obtain converged solutions regardless 
of boundary conditions. For A= 0° (see fig. 3(a), direction of maximum flexural stiff- 
ness alined with flow), converged flutter values are obtained for only four terms in the 
D 1 direction for all values of for which calculations were made. For 
the A*' ^ values converge to the exact value obtained from reference 17. The converged 
\l/3 value for = 10 is slightly higher than the exact value shown for free edge con- 
ditions = oj (see ref. 7). For A= 90° (see fig. 3(b)), the situation is quite different, 
as converged results are seen to be strongly dependent on the boundary conditions. Con- 
verged flutter values are obtained with four terms in the D 2 direction for = 10, with 
eight terms in the D 2 direction for K d = 100, and with 20 terms in the D 2 direction 
for Kjj = 500. For K<j = 00 (simply supported edge), the flutter solution is not converged 
when even 60 terms are used; xV 2 is considerably below the exact value and, further, 
convergence is extremely slow. Thus, for orthotropic panels with simply supported edges 
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with flow normal to the direction of maximum flexural stiffness, a large number of modal 
terms are needed to obtain converged solutions. Similar trends have also been found to 
be true for panels with clamped edges (see ref. 13). As K d decreases, however, the 
number of terms required for convergence decreases rapidly. When using only the first 
two terms in the D2 direction (n = 1,2) for K d = 100 and 500 (see fig. 3(b)), X 1 / 3 is 
seen to be considerably higher than the converged values. For these values of K d , 
higher modes coalesce for the lowest value of X 1 / 3 , and the absence of these modes in 
the two-term solution results in erroneous values. 

Convergence results for flow at A= 10° are shown in figure 3(c) where X 1 / 3 is 
plotted against the number of modal terms in the D2 direction. Flutter points are shown 
for Dj modal terms from 1 to 5. For A = 10°, one term in the Dj direction still 
gives converged solutions for K d = 10, 100, and 500 and for 4, 8, and 16 terms in the 
D2 direction, respectively. For K d = 00 (simply supported), however, one term in the 
Dj direction results in a diverging solution which becomes progressively worse as more 
terms in the D 2 direction are included in the calculations. When three terms are used 
in the Dj direction, only four terms are required in the D 2 direction to give converged 
flutter values for K d = °®, 

The variation in the modal terms required for convergence is shown in figure 4 
where X*/ 3 is plotted as a function of flow angle for Djg/Dj = 0.5, Dj 2 / D 2 “ and 
(a/b) 2 = 1.0. The curves are the flutter boundaries for the values of K d shown. The 
numbers in parentheses indicate the number of terms in the Dj and D2 directions 
required to obtain convergence within ±2 percent. The K d = 00 curve is not shown for 
A > 10° owing to the large amount of computer time required for convergence. The data 
show that, even for small flow angles, considerably more D2 terms are required for 
convergence than Dj terms. For A= 30° only one term in the Dj direction is 
required for all K d values shown, and the number of terms required in the D2 direction 
for A = 30° is adequate for all higher values of A. 

Flutter Frequencies and Mode Shapes 

Further insight into the modal terms required for convergence can be gained 
through examination of the influence of flow angle on frequencies and mode shapes. 

Panel frequencies for Di 2 / D i “ 0-5, d 12/^2 85 500, an< * (a/b) 3 = 1-0 with Kd = 100 
are shown in figure 5 as a function of X for flow angles of 0°, 2°, 15°, and 90°. 

For A = 0° (fig. 5(a)), increasing X results in the coalescence of the frequencies of 
modes m = 1 and 2 and n = 1 at X cr =118. Further increases in X result in a pair 
of frequencies becoming complex conjugates as shown by the dashed curve. Although the 
(1,1) and (1,2) frequencies are identical at X =75, panel instability does not occur. For 
A= 2° (fig. 5(b)), however, coalescence of the (1,1) and (1,2) modes does occur at 
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Xc r = 70 and remains unstable for further increases in X. In reference 10, where two 
terms were used in each direction, a stable region was found at values of X greater 
than the initial instability value. Such a stable region clearly does not exist if a sufficient 
number of modes are used in the analysis. These data also suggest that, for small 
increases in A from zero, large reductions in critical dynamic pressure occur. This 
behavior must be taken into account in the design of wavy (highly orthotropic) surface 
panels since some cross flow will be experienced in actual flight. 

For A= 15° (fig. 5(c)), the (1,3) and (1,4) modes coalesce at a X cr value of 
approximately 17, which represents a large drop from that obtained at A - 2°. For 
A= 90° (fig. 5(d)), the (1,3) and (1,4) modes still coalesce but for values of =5. 

The flutter mode shapes for = 0*5, Dj 2/1*2 = 500, (a/b)^ = 1, and 

= 100 are shown in figure 6 for flow angles of 0°, 2°, 15°, and 90°. For 
A = 0° (fig. 6(a)), the mode shape has two half-waves in the direction of maximum 
panel flexural stiffness and one half-wave in the direction of minimum panel flexural 
stiffness. The maximum deflection occurs on the downstream portion of the panel. At 
A = 2° (fig. 6(b)), the flutter mode shape is considerably different from that at A = 0°, 
with the maximum deflection shifting in the direction of minimum flexural stiffness. At 
a flow angle of 15° (fig. 6(c)) the flutter mode shape is complicated but basically is com- 
posed of three half-waves in the direction of minimum flexural stiffness and one half-wave 
in the opposite direction. Further increases in A have little additional effect on this 
mode shape (see fig. 6(d) for A = 90°). 


Parametric Studies 

Effects of flow angle .- The effects of flow angle on flutter are shown in figure 7 
where flutter boundaries are presented as a function of flow angle for D 12/ D l “ °* 05 
and (a/b)^ = 10. Boundaries are shown for = 10 (fig. 7(a)), K d = 100 (fig. 7(b)), 
and K(j = 00 (fig. 7(c)). The curves of each figure are for values of of 2, 20, 

50, and 500. The actual calculations are indicated by the symbols on the curves. For 
K(j = 10 (fig. 7(a)) and Dj 2/^2 = 500, an order-of-magnitude reduction in the flut- 
ter dynamic pressure occurs over the first 10°, with less severe reductions there- 
after. For smaller values of Dj2/ D 2 (corresponding to increases in D2^, the reduc- 
tion in flutter dynamic pressure with flow angle is less and is spread out over a wider 
flow-angle range. For £>12/^2 = 2, even an increase in the flutter dynamic pressure is 
obtained with increasing flow angles over most of the flow-angle range. This latter 
behavior is similar to that expected for an isotropic panel with a similar a/b ratio. 
Two points are shown at A= 0° for 2/^2 = ^00 and are connected by a dashed line. 
The upper point is the critical flutter value at A= 0° and is the same for all values 
of 2/1*2* The lower point is the value of X where the two lowest frequencies cross 
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(see fig. (5)). Increasing the flow angle from zero results in the critical flutter X 
dropping from the higher value at A = 0° to the lower curve for a very small A (for 
example, A = 1°). Thus, for all practical purposes, the lower values should be used, 
even for A = 0°. 

Increasing the deflection spring stiffness, K d , to 100 (fig. 7(b)) and 00 (fig. 7(c)) 
results in flutter boundaries similar to those obtained for K d = 10 but higher throughout 
the flow- angle range. 

Flutter boundaries for = 0.25 and 0.5 are shown in figures 8 and 9, 

respectively, and have trends similar to those discussed previously for D^/^l = 0-05. 
Panels with large values 3X6 more stable, however, as indicated by higher 

values of X throughout the flow-angle range. For Dj2/ D l = 0.25 ( fi S‘ 311(1 A= 8 °» 
note that X(. r is lower for K d = 100 than for K d = 10. This unusual trend has been 
shown previously by an exact analysis in reference 7 and is attributed to a change in the 
flutter mode. With K d = 10 the lower vibration modes approach those of rigid-body 
translation and rotation and the panel becomes more resistant to flutter. At D 12 / / D l = 0,5 
(fig. 9), this phenomenon does not occur. 

Effects of flexural stiffness .- The data presented in figures 7 to 9 are replotted 
in figure 10 as a function of D 1 2/D 1 for flow angles of 0°, 5°, 10°, and 20°. Flut- 
ter values nondimensionalized by X cr . (the value of X for a simply supported 

(Kd= w ) _ 

panel at the flow angle in question) are shown for (a/b) 2 = 10, K d = 100, and D j2/ D 2 
of 2, 20, 50, and 500. In figure 10(a), the boundaries for A= 0° are independent of 
values of Dj 2 /t>2 and are, therefore, functions only of K d . Increasing values of stiff- 
ness ratio, ^ 12 /^V resu5 - t in flutter boundaries approaching those for the simply sup- 
ported panel (K d = =°). The flutter boundary is also shown for K d = 10 because of its 
peculiar trend. The unusual behavior for K d = 10 is due to the phenomenon that was 
noted previously in the discussion of figure 8 for D^/^l = 0.25. Flutter values are 
obtained for K d = 10 that are even higher than for simply supported panels for the range 
of ^>12/^1 f rom 0*1 to 0.2. This behavior is consistent with results obtained from an 
exact analysis presented in reference 7. 

Flutter boundaries for A> 0° are seen to be dependent on both 311,1 

D 12A > 2* ^ or smaJ1 - values of £>12/^1 (* or exam Pl e » Di 2 /Dj = 0.05 j, increases in D 2 
2/^2 decreases) have little effect on X; however, for large values of D 1 2/° 1 
^for example, = 0.5) increases in Dg have a large effect on X. Note also that 

the influence of is strongly dependent on the flow angle; for example, with 

A = 10° and I>i2^> 2 > 50, changes in Dj have negligible inf luence on flutter. Thus, 
the curves tend to identify the range of influence of stiffnesses and D 2 and flow 
angle on panel flutter. 
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Effects of length-width ratio .- Flutter boundaries for = 0-5 and 

(a/b)^ = 1.0 are presented in figure 11. The effects of panel length-width ratio on the 
flutter boundaries may be seen by comparing the boundaries given in figure 11 with those 
presented in figure 9. Reducing the panel length-width ratio has a large destabilizing 
effect as indicated by considerably lower values of X throughout the flow-angle and 
K,j range. The destabilizing effect is more pronounced for elastically supported panels 
than for simply supported panels. 

Effects of Aerodynamic and Structural Damping 

The effects of aerodynamic and structural damping on the flutter of unstressed 
orthotropic panels are shown in figure 12. Calculations are made for ~ 

D 12/ D 2 = 500 > (a/* 5 ) 2 = 1«°> = 100, and for values of g a = 0.2 and g b = 0.05. The 

membrane structural damping coefficient, g m , is not considered since it has an effect 
only if there are in-plane forces acting on the panel. As can be seen, aerodynamic 
damping has a significant stabilizing effect at all flow angles except zero, whereas struc- 
tural damping has negligible effect. 

CONCLUDING REMARKS 

A theoretical panel flutter analysis and computer program was developed which is 
capable of analyzing an orthotropic panel with in-plane loads, at various angles of cross 
flow, and with two opposite edges mounted on flexible supports and the other two edges 
simply supported. The resulting modal analysis uses linear piston-theory aerodynamics 
and includes both aerodynamic and structural damping. The number of modal terms 
used in each direction to analyze the panel can be varied to insure that converged solu- 
tions are obtained. 

The inclusion of cross-flow angle as a variable in a panel flutter analysis requires 
that the panel geometry/ stiffness properties be defined. Thus, flutter calculations must 
be made for each specific panel. Calculations made for typical unstressed orthotropic 
panels with deflectional spring supports indicate that fewer modal terms are required for 
converged solutions for elastically supported panels than for simply supported panels. 

The calculations further show that an order-of-magnitude reduction in flutter 
dynamic pressure is possible with only small changes in flow angle. Also, the effects of 
flow angle are more pronounced for elastically supported panels than for simply sup- 
ported panels. In addition, the results indicate that the influence of flexural stiffness, 

Dj, on flutter of elastically supported panels diminishes with small increases in flow 
angle from zero and that the range of influence of flexural stiffness can be identified as 
a function of flow angle. 
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Figure 1.- Typical metallic thermal protection system 


L-73-5498.1 
















28 







T~-\ 




Figure 6.- Concluded. 
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Figure 10.- Influence of panel flexural stiffness on flutter, (a/b) = 10, g a =0, g b = 0, and g 






Figure 12.- Effects of aerodynamic and structural damping on flutter. 
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